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Abstract 

Within the framework of a dinuclear system (DNS) model, the evaporation-residue excitation 
, ^/ functions and the quasi-fission mass yields in the 48 Ca induced fusion reactions are investigated 
CN) ■ systematically and compared with available experimental data. Maximal production cross sec- 

> : 

tions of superheavy nuclei based on stable actinide targets are obtained. Isotopic trends in the 
production of the superheavy elements Z=110, 112-118 based on the actinide isotopic targets 
CT) ■ are analyzed systematically. Optimal evaporation channels and combinations as well as the cor- 

o : 

responding excitation energies are proposed. The possible factors that influencing the isotopic 

O ■ 

dependence of the production cross sections are analyzed. The formation of the superheavy nuclei 

^ ■ based on the isotopes U with different projectiles are also investigated and calculated. 
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1 Introduction 

The synthesis of heavy or superheavy nuclei is a very important subject in nuclear physics mo- 
tivated with respect to the island of stability which is predicted theoretically, and has obtained 
much experimental research with the fusion-evaporation reactions [H [2]. The existence of the 
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superheavy nucleus (SHN) {Z > 106) is due to strong binding shell effects against the large 
Coulomb repulsion. However, the shell effects get reduced with increasing the excitation energy 
of the formed compound nucleus. Combinations with a doubly magic nucleus or nearly magic 
nucleus are usually chosen owing to the larger reaction Q values. Reactions with 208 Pb or 209 Bi 
targets were first proposed by Oganessian et al. to synthesize SHN [3J. Six new elements with 
Z=107-112 were synthesized in cold fusion reactions for the first time and investigated at GSI 
(Darmstadt, Germany) with the heavy-ion accelerator UNILAC and the SHIP separator [fl H]. 
Recently, experiments on the synthesis of element 113 in the 70 Zn+ 209 Bi reaction have been per- 
formed successfully at RIKEN (Tokyo, Japan) [5J. However, it is difficulty to produce heavier 
SHN in the cold fusion reactions because of the smaller production cross sections that are lower 
than 1 pb for Z > 113. Other possible ways to produce SHN are very needed to be investigated in 
experimentally and theoretically. Recently, the superheavy elements Z=113-116, 118 were synthe- 
sized at FLNR in Dubna (Russia) with the double magic nucleus 48 Ca bombarding actinide nuclei 
[HI El IB]- New heavy isotopes 259 Db and 265 Bh have also been synthesized at HIRFL in Lanzhou 
(China) [9]. Further experimental works are necessary in order to testify the new synthesized 
SHN. A reasonable understanding of the formation of SHN in the massive fusion reactions is still 
a challenge for theory. 

In accordance with the evolution of two heavy colliding nuclei, the dynamical process of the 
compound nucleus formation and decay is usually divided into three reaction stages, namely the 
capture process of the colliding system to overcome the Coulomb barrier, the formation of the 
compound nucleus to pass over the inner fusion barrier, and the de-excitation of the excited 
compound nucleus by neutron emission against fission. The transmission in the capture process 
depends on the incident energy and relative angular momentum of the colliding nuclei, which is 
the same as that in the fusion of light and medium mass systems. The complete fusion of the 
heavy system after capture in competition with quasi-fission is very important in the estimation 
of the SHN production. The concept of the "extra-push" energy explains for the fusion of two 
heavy colliding nuclei in the macroscopic dynamical model [TQl CEE]- At present it is still difficult 
to make an accurate description of the fusion dynamics. After the capture and the subsequent 
evolution to form the compound nucleus, the thermal compound nucleus will decay by the emission 
of light particles and 7 rays against fission. The three stages will affect the formation of evapo- 
ration residues observed in laboratories. The evolution of the whole process of massive heavy-ion 
collisions is very complicated at near-barrier energies. Most of the theoretical methods on the 
formation of SHN have a similar viewpoint in the description of the capture and the de-excitation 
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stages, but there are different description of the compound nucleus formation process. There are 
mainly two sorts of models, depending on whether the compound nucleus is formed along the 
radial variable (internuclear distance) or by nucleon transfer in a touching configuration which is 
usually the minimum position of the interaction potential after capture of the colliding system. 
Several transport models have been established to understand the fusion mechanism of two heavy 
colliding nuclei leading to SHN formation, such as the macroscopic dynamical model [101 |TT], the 
fluctuation-dissipation model [12], the concept of nucleon collectivization [13] and the dinuclear 
system model [Tl"l ITS] . Recently, the improved isospin-dependent quantum molecular dynamics 
(ImlQMD) model was also proposed to investigate the fusion dynamics of SHN [T6J [17]. With 
these models experimental data can be reproduced to a certain extent, and some new results have 
been predicted. However, these models differ from each other, and sometimes different physical 
ideas are used. 

Further improvements of these models have to be made. Here we use a dinuclear system (DNS) 
model [151 HE], m which the nucleon transfer is coupled with the relative motion by solving a set 
of microscopically derived master equations, and a barrier distribution of the colliding system 
is introduced in the model. We present a new and extended investigation of the production of 
superheavy nuclei in the 48 Ca induced fusion reactions and in other combinations. 

In Section 2 we give a simple description on the DNS model. Calculated results of fusion 
dynamics and SHN production are given in Section 3. In Section 4 conclusions are discussed. 

2 Dinuclear system model 

The dinuclear system [19] is a molecular configuration of two touching nuclei which keep their 
own individuality [2] . Such a system has an evolution along two main degrees of freedom: (i) the 
relative motion of the nuclei in the interaction potential to form the DNS and the decay of the 
DNS (quasi-fission process) along the R degree of freedom (internuclear motion), (ii) the transfer 
of nucleons in the mass asymmetry coordinate 7] = (Ai — A 2 )/(Ai + A 2 ) between two nuclei, 
which is a diffusion process of the excited systems leading to the compound nucleus formation. 
Off-diagonal diffusion in the surface (Ai, R) is not considered since we assume the DNS is formed 
at the minimum position of the interaction potential of two colliding nuclei. In this concept, the 
evaporation residue cross section is expressed as a sum over partial waves with angular momentum 
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°E R {E c . m ) = T-= E ( 2J + l)T(E c . m ., J)P C N(E c . m ., J)W sur (E c . m ., J). (1) 
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Here, T(E c _ m _, J) is the transmission probability of the two colliding nuclei overcoming the Coulomb 
potential barrier in the entrance channel to form the DNS. In the same manner as in the nucleon 
collectivization model [13], the transmission probability T is calculated by using the empirical cou- 
pled channel model, which can reproduce very well available experimental capture cross sections 
[33J [15]. The Pcn is the probability that the system will evolve from a touching configuration 
into the compound nucleus in competition with quasi-fission of the DNS and fission of the heavy 
fragment. The last term is the survival probability of the formed compound nucleus, which can be 
estimated with the statistical evaporation model by considering the competition between neutron 
evaporation and fission [15J. We take the maximal angular momentum as J max = 30 since the 
fission barrier of the heavy nucleus disappears at high spin [20] . 

In order to describe the fusion dynamics as a diffusion process in mass asymmetry, the analyti- 
cal solution of the Fokker-Planck equation [H] and the numerical solution of the master equations 
[2T| 122] have been used, which were also used to treat deep inelastic heavy- ion collisions [23] . 
Here, the fusion probability is obtained by solving a set of master equations numerically in the 
potential energy surface of the DNS. The time evolution of the distribution function P(A\, Ei,t) 
for fragment 1 with mass number A\ and excitation energy E\ is described by the following master 
equations (HI |2T] , 

dP{A x E x ,t) = j2w AuA[ (t) [d Al P(A' l7 E[ : t) - dj,P{A x ,E 1: t) - 
at A , 

A 9/ (e(t)) + A fts (Q(t))} P(A 1 ,E 1 ,t). (2) 

Here HOii.A'j is the mean transition probability from the channel {A\,E\) to (A^E^), and d Al 
denotes the microscopic dimension corresponding to the macroscopic state (A\,E\). The sum is 
taken over all possible mass numbers that fragment A[ may take (from to A = A\ + A 2 ), but 
only one nucleon transfer is considered in the model with A\ = A\±l. The excitation energy E\ is 
the local excitation energy e\ with respect to fragment A\, which is determined by the dissipation 
energy from the relative motion and the potential energy of the corresponding DNS and will be 
shown later in Eqs.(8) and (9). The dissipation energy is described by the parametrization method 
of the classical deflection function [2~31 [25] . The motion of nucleons in the interacting potential is 
governed by the single-particle Hamiltonian [T5J [21]: 

H(t) = H (t) + V{t) (3) 



with 



H (t) = EE^WiW^W. 
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v(t) = E E = E v K , K ,(t). (4) 

K,K'a K ,f3 K , K,K' 

Here the indices K,K' (K,K' = 1,2) denote the fragments 1 and 2. The quantities e UK and 
Ua K ,(3 K , represent the single particle energies and the interaction matrix elements, respectively. 
The single particle states are defined with respect to the centers of the interacting nuclei and are 
assumed to be orthogonalized in the overlap region. So the annihilation and creation operators 
are dependent on time. The single particle matrix elements are parameterized by 

1 ( £ajc(*) - e P K X t )\ 



*W,/3 K ,(*) = U KiK ,{t) I exp 
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which contain some parameters UK,K'{t) and A^^/(t). The detailed calculation of these parame- 
ters and the mean transition probabilities were described in Refs. [T5| |2T] . 

The evolution of the DNS along the variable R leads to the quasi-fission of the DNS. The 
quasi-fission rate A q f can be estimated with the one-dimensional Kramers formula [261 127] : 

Here the quasi-fission barrier is counted from the depth of the pocket of the interaction potential. 



The local temperature is given by the Fermi-gas expression = ye* /a corresponding to the local 
excitation energy e* and level density parameter a = A/12 MeV~ l . In Eq.(6) the frequency uj B if 
is the frequency of the inverted harmonic oscillator approximating the interaction potential of two 
nuclei in R around the top of the quasi-fission barrier, and uj is the frequency of the harmonic 
oscillator approximating the potential in R around the bottom of the pocket. The quantity T, 
which denotes the double average width of the contributing single-particle states, determines the 
friction coefficients: jw = \nu>, with /i^/ being the inertia tensor. Here we use constant values 
T = 2.8 MeV, hu Bqf = 2.0 MeV and Tiuo = 3.0 MeV for the following reactions. The Kramers 
formula is derived with the quasi- stationary condition of the temperature 0(t) < B q j(Ai, A2). 
However, the numerical calculation in Ref . [27] indicated that Eq. (6) is also useful for the condition 
of Q(t) > B q f(Ai, A 2 ). In the reactions of synthesizing SHN, there is the possibility of the fission 
of the heavy fragment in the DNS. Because the fissility increases with the charge number of 
the nucleus, the fission of the heavy fragment can affect the quasi-fission and fusion when the 
DNS evolves towards larger mass asymmetry. The fission rate A^ iS can also be treated with the 



one-dimensional Kramers formula 



where the u g , s , and cj/ are the frequencies of the oscillators approximating the fission-path potential 
at the ground state and on the top of the fission barrier for nucleus A\ or A 2 (larger fragment), 
respectively. Here, we take Tnrf g . s . = huf = 1.0 MeV, T = 2 MeV. The fission barrier is calculated 
as the sum of a macroscopic part and the shell correction energy used in Refs. [151 I2H]- The fission 
of the heavy fragment does not favor the diffusion of the system to a light fragment distribution. 
Therefore, it leads to a slight decrease of the fusion probability. 

In the relaxation process of the relative motion, the DNS will be excited by the dissipation of the 
relative kinetic energy. The excited system opens a valence space Aek in fragment K(K = 1, 2), 
which has a symmetrical distribution around the Fermi surface. Only the particles in the states 
within this valence space are actively involved in excitation and transfer. The averages on these 
quantities are performed in the valence space: 
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where the e* is the local excitation energy of the DNS, which provides the excitation energy for 
the mean transition probability. There are Nk = QkAek valence states and mx = Nk/2 valence 

I N x \ I N 2 \ 

nucleons in the valence space Aek, which gives the dimension d(m\, m 2 ) — I j • The 

y mi J y m 2 J 

local excitation energy is defined as 

e* = E x - (U(A 1 , A 2 ) - U(A P , A T )) . (9) 

Here the U(A 1 ,A 2 ) and U(Ap, A?) are the driving potentials of fragments Ai, A 2 and fragments 
Ap, At (at the entrance point of the DNS), respectively. The detailed calculation of the driving 
potentials can be seen in Ref. [T5]. The excitation energy E x of the composite system is converted 
from the relative kinetic energy loss, which is related to the Coulomb barrier B [2H] and determined 
for each initial relative angular momentum J by the parametrization method of the classical 
deflection function [2^ [25] . So E x is coupled with the relative angular momentum. 

After reaching the reaction time in the evolution of P(Ai, Ei, t), all those components on 
the left side of the B.G. (Businaro-Gallone) point contribute to the formation of the compound 
nucleus. The hindrance in the diffusion process by nucleon transfer to form the compound nucleus 
is the inner fusion barrier Bf us , which is defined as the difference of the driving potential at the 



B.G. point and at the entrance position. Nucleon transfers to more symmetric fragments undergo 
quasi-fission. The formation probability of the compound nucleus at the Coulomb barrier B (here 
a barrier distribution f(B) is considered) and angular momentum J is given by 

A B g 

P C N(E c . m .,J,B) = J2 P{A 1: E l: r int {E c _ m . : J : B)). (10) 



Here the interaction time Ti nt (E c , m , } J, B) is obtained using the deflection function method [30] . 
which means the time duration for nucleon transfer from the capture stage to the formation of 
the complete fused system with the order of 10 -20 s. We obtain the fusion probability as 

Pc N (E c . m ., J)= J f(B)P CN (E c . m „ J, B)dB, (11) 



where the barrier distribution function is taken in asymmetric Gaussian form [13| 15J. So the 
fusion cross section is written as 
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cT f us(E c . m .) = — — ]T(2J + l)T(E c . m ., J)P CN (E c . m ., J). (12) 
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The survival probability of the excited compound nucleus cooled by the neutron evaporation 
in competition with fission is expressed as follows: 

WM,*,J) = «»,an ( r^.^ + r'U.j) ),- < 13 > 

where the Eq N , J are the excitation energy and the spin of the compound nucleus, respectively. 
The E* is the excitation energy before evaporating the ith neutron, which has the relation 

E* +1 = E* - B? - 2T h (14) 

with the initial condition E^ = Eq N . The energy B™ is the separation energy of the ith neu- 
tron. The nuclear temperature T; is given by E* = aTf — Ti with the level density parameter a. 
P(Eq N ,x, J) is the realization probability of emitting x neutrons. The widths of neutron evapo- 
ration and fission are calculated using the statistical model. The details can be found in Ref. [TB] . 
The level density is expressed by the back-shifted Bethe formula [31J with the spin cut-off model 
as 

2^ + 1 -1/4,™ _ AV -5/4 r „„ro Lr^ a^ , ( J + V2) 2 



p(E\ J) = K^K^—^-a-^iE* - A)- 5 / 4 exp[2 Ja(E* - A)] exp[- v ' ' ' ], (15) 
24 v la 6 lo 

where the K rot and K vi i } are the coefficients of the rotational and vibrational enhancements. The 
pairing energy is given by 
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in MeV(x=-l, and 1 for odd-odd, odd-even and even-even nuclei, respectively). The spin cut-off 
parameter is calculated by the formula: 

a 2 = T( r . b /h\ (17) 

where the rigid-body moment of inertia has the relation £ r .& = OAMR 2 with the mass M and 
the radius R of the nucleus. The level density parameter is related to the shell correction energy 
E sh (Z, N) and the excitation energy E* of the nucleus as 

a{E\ Z, N) = a(A)[l + E sh (Z, N)f(E* - A)/(E* - A)]. (18) 

Here, a(A) = aA + (3A 2 ^ 3 b s is the asymptotic Fermi-gas value of the level density parameter at 
high excitation energy. The shell damping factor is given by 

/(£*) = l-exp(- T £*) (19) 

with j — a/ (eA 4 / 3 ). All the used parameters are listed in Table 1. In Fig.l we give the level density 
parameters of different nuclides at the ground state calculated by using Eq.(18) and compared 
them with two empirical formulas a(A) = A/8, and A/12. It can be seen that the strong shell 
effects appear in the level density. 

With this procedure introduced above, we calculated the angular momentum dependence of 
the capture, fusion and survival probabilities as shown in Fig. 2 for the reaction 48 Ca+ 208 Pb at 
incident energies 172.36 MeV and 192.36 MeV, respectively. The values of the three stages decrease 
obviously with increasing the relative angular momentum. So in the following estimation of the 
production cross sections, we cut off the maximal angular momentum at J ma x — 30, which is taken 
as the same value that used in the cold fusion reactions [18J. 

3 Results and discussions 

3.1 Fusion-fission reactions and quasi-fission mass yields 

As a test of the parameters for the estimation of the transmission of two colliding nuclei and the 
de-excitation of the thermal compound nucleus, we analyzed the fusion-fission reactions for the 
selected systems shown in Fig. 3 assuming P C n = 1- The capture and evaporation residue cross 
sections are compared with the available experimental data [221 E21 EH 135] . For these systems 
the quasi-fission does not dominate in the sub-barrier region, which also means that Pcn ~ 1- 



The evaporation residues are mainly determined through the capture of the light projectile by the 
target nucleus and the survival probabilities of the formed compound nucleus. The experimental 
data can be reproduced rather well within the error bars. Some discrepancies may come from 
the quasi-fission in the above barrier region and from the input quantities, such as the neutron 
separation energy, shell correction and mass. The rotational and the vibrational enhancement in 
the level density can also affect the survival probabilities of the excited compound nucleus [36] . 
Here we take unity for both coefficients as shown in Table 1 because the height of the fission barrier 
is also sensitive to the survival of the compound nucleus by fitting the experimental evaporation 
residue excitation functions in the fusion-fission reactions. 

Since the electrostatic energy of the composite systems formed by two heavy colliding nuclei 
is very large, so although the two nuclei may be captured by the nuclear potential, they almost 
always separate after mass transfer from the heavier nucleus to the lighter one rather fusing. This 
process is called quasi-fission [33, [38], which is the main feature in the massive fusion reactions 
and can inhibit fusion by several degree of freedom. Recently, experiment has performed nice 
works by measuring the quasi-fission and fusion-fission mass yields [52]. In the DNS model, the 
quasi-fission mass yields are expressed as [26J 




In Fig. 4 we show a comparison of the calculated quasi-fission mass yields and the experimental data 
for the two 48 Ca induced reaction systems. The trends of the distribution can be reproduced by 
the DNS model. At the domain of the medium-mass fragments Ai=Acw/2-30~Acw/2-|-30, The 
experimental data are higher than the calculated values, which may be come from the contribution 
of the fusion-fission fragments. 

3.2 Evaporation residue cross sections 

The evaporation residues observed in laboratories by the consecutive a decay are mainly produced 
by the complete fusion reactions, in which the fusion dynamics and the structure properties of the 
compound nucleus affect their production. Within the framework of the DNS model, we calculated 
the evaporation residue cross sections producing SHN Z=110, 112, 113, 115 with 232 Th, 238 U, 237 Np 
and 243 Am targets in the 48 Ca induced reactions as shown in Fig.5, and compared them with the 
Dubna data [7J HUJ UT] as well as with the recent GSI data [12] for 238 U targets in the 3n channel. 
Compared with the Dubna data for the system 48 Ca+ 238 U, the GSI results show that the formation 
cross sections in the 3n channel have a slight decrease at the same excitation energy, which is in 
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a good agreement with our calculated results. The calculations were carried out before getting 
the experimental data [JT] for the reaction 48 Ca+ 237 Np, and a good agreement with the data is 
also found [33]. The excitation energy of the compound nucleus is obtained by Eq N = E c _ m , + Q, 
where the E c _ m _ is the incident energy in the center-of-mass system. The Q value is given by 
Q = AMp + AMt — AMq, and the corresponding mass excesses AMj (i = P, T, C) are taken the 
data from Ref. [44J for the projectile, target and compound nucleus denoted with the symbols 
P, T and C, respectively. Usually, the neutron-rich projectile-target combinations are in favor of 
synthesizing SHN experimentally, which can enhance the survival probability W sur in Eq.(l) of 
the formed compound nucleus because of the smaller neutron separation energy. Differently to 
the cold fusion reactions |18j . the maximal production cross sections from Ds to 115 especially 
in the 2n-5n channels are not changed much although the heavier SHNs are synthesized. Within 
the error bars the experimental data can be reproduced rather well. With the same procedure, 
we analyzed the evaporation residue excitation functions with targets 242 ' 244 Pu and 245,248 Cm that 
are used to synthesize the superheavy elements Z=114 and 116 in Dubna [1Q1 US] (Fig. 6). Our 
calculations show that the target 244 Pu has a larger production cross section than 242 Pu because 
of the larger survival probability. In Fig. 7 we also calculated the evaporation residue excitation 
functions to synthesize superheavy elements Z=117-120 using the actinide isotopes with longer 
half-lives 247 Bk, 249 Cf, 254 Es and 257 Fm. The 3n evaporation channel with an excitation energy 
of the formed compound nucleus around 30 MeV is favorable to produce SHN with Z>117 by 
using the actinide targets. Within the error bars, the positions of the maximal production cross 
sections are in good agreement with the available experimental results. Similar calculation of 
the evaporation residue excitation functions was also reported in Ref. [IB]. The spectrum form 
of evaporating neutrons is mainly determined by the survival probability, in which the neutron 
separation energy and the shell correction play a very important role in the determination of the 
value. We considered the angular momentum influence in the calculation of the level density, but 
did not include it in the estimation of the fission barrier of the thermal compound nucleus. As 
pointed out in section 1, the fission barrier of SHN decreases rapidly with increasing excitation 
energy of the compound nucleus, where the rotation of the system affects the height of the barrier 
and also influences other crucial quantities such as the level density etc. 

In Fig. 8 we show a comparison of the calculated maximal production cross sections of super- 
heavy elements Z=102-120 in the cold fusion reactions by evaporating one neutron, in the 48 Ca 
induced reactions with actinide targets by evaporating three neutrons, and the experimental data 
[U [2], HI H?]. The production cross sections decrease rapidly with increasing the charge number of 
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the synthesized compound nucleus in the cold fusion reactions, such as from 0.2 \ib for the reaction 
48 Ca+ 208 Pb to 1 pb for 70 Zn+ 208 Pb, and even below 0.1 pb for synthesizing Z>113 [18j. It seems 
to be difficult to synthesize superheavy elements Z>113 in the cold fusion reactions at the present 
facilities. The calculated results show that the 48 Ca induced reactions have smaller production 
cross sections with 232 Th target, but are in favor of synthesizing heavier SHN (Z>113) because of 
the larger cross sections. The experimental data also give such trends. In the DNS concept, the 
inner fusion barrier increases with reducing mass asymmetry in the cold fusion reactions, which 
leads to a decrease of the formation probability of the compound nucleus. However, the 48 Ca 
induced reactions have not such increase of the inner fusion barrier for synthesizing heavier SHN. 
Because of the larger transmission and the higher fusion probability, we obtain larger production 
cross sections for synthesizing SHN (Z>113) in the 48 Ca induced reactions although these reactions 
have the smaller survival probability than those in the cold fusion reactions. It is still a good way 
to synthesize heavier SHN by using the 48 Ca induced reactions. Of course, further experimental 
data are anticipated to be obtained in the future. However, the actinide targets are difficulty to 
be handled in experiments synthesizing heavier SHN. 

3.3 Isotopic dependence of the production cross sections 

Recent experimental data show that the production cross sections of the SHN depend on the 
isotopic combination of the target and projectile in the 48 Ca induced fusion reactions. For example, 
the maximal cross section in the 3n channel is 3.7± 3 /g pb for the reaction 48 Ca+ 245 Cm at the 
excitation energy 37.9 MeV; however, it is 1.2 pb for the reaction 48 Ca+ 248 Cm although the later 
is a neutron-rich target |E, 40J. The isotopic trends of the production cross sections were also 
observed and investigated in cold fusion reactions [HI [18]. Further investigations on the isotopic 
trends in the 48 Ca induced reactions are very necessary for predicting the optimal combinations, 
excitation energies (incident energies) and evaporation channels in the synthesis of SHN. In Fig. 9 
we show the calculated isotopic trends in producing superheavy elements Z=110, 112 with the 
isotopic actinides Th and U in the 3n channels, and compare them with the available experimental 
data performed in Dubna [10] (squares with error bars) and at GSI [12] (circles with error bars). 
The results show that the targets 230 Th in the 4n channel and 235 > 238 U in the 3n channel have the 
largest cross sections. The isotopic trends in synthesizing Z=113-116 with the actinide targets 
Np, Pu, Am and Cm are also calculated systematically, and compared with the existing data 
measured in Dubna [7J HUJ H5] and the results of Adamian et al. [H] for the Pu isotopes as shown 
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in Fig.10 and Fig.ll. The isotopes 237 Np, 241 Pu, 242,243 Am and 245 ' 247 Cm in the 3n channels, and 
244 Pu in the 4n channel as well as the isotope 250 Cm are suitable for synthesizing SHN. Except 
for the 244 Pu, our calculated cross sections are smaller than the ones of the Adamian et al. In 
the DNS model, the isotopic dependence of the production cross sections is mainly determined by 
both the fusion and survival probabilities. Of course, the transmission probability of two colliding 
nuclei can also be affected since the isotopes have initial quadrupole deformations. With the 
same procedure, we analyzed the dependence of the production cross sections on the isotopes Bk 
and Cf in the 3n channels for synthesizing the superheavy elements Z=117, 118 and compared 
them with the available experimental data [8] shown in Fig. 12. The results show that the targets 



248 ' Bk and 251 ' 252 Cf are favorable for synthesizing the superheavy elements Z=117 and 118. The 
corresponding excitation energies are also given in the figures. 

In Fig. 13 we show the dependence of the inner fusion barrier, the fission barrier of the com- 
pound nucleus, and the neutron separation energies of evaporating 3n and 4n on the mass numbers 
of the isotopic targets Cm in the 48 Ca induced reactions. It is obvious that the combinations with 
the isotopes 245,247 Cm have smaller inner fusion barriers, higher fission barriers and smaller 3n 
separation energies, which result in larger production cross sections producing the superheavy 
element Z=116. Although the lower fission barrier for the isotope 250 Cm, it gives the smaller 
inner fusion barrier and neutron separation energies, which also leads to the larger cross sections 
in the 3n and 4n channels as shown in Fig.ll. The shell correction and the neutron separation 
energies are taken from Ref. [H]. When the neutron number of the target increases, the DNS 
gets more asymmetrical and the fusion probability increases if the DNS does not consist of more 
stable nuclei (such as magic nuclei) because of a smaller inner fusion barrier. A smaller neutron 
separation energy and a larger shell correction lead to a larger survival probability. The compound 
nucleus with closed neutron shells has a larger shell correction energy and a larger neutron sep- 
aration energy. The neutron-rich actinide target has larger fusion and survival probabilities due 
to the larger asymmetric initial combinations and smaller neutron separation energies. But such 
actinide isotopes are usually unstable with smaller half-lives. With the establishment of the high 
intensity radioactive-beam facilities, the neutron-rich SHN may be synthesized experimentally, 
which approaches the island of stability. 
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3.4 238 U based reactions 

The uranium is the heaviest element existing in the nature. It has a larger mass asymmetry 
constructed as a target in the fusion reactions with the various neutron-rich light projectiles. The 
isotope 238 U is the neutron-richest nucleus in the U isotopes and often chosen as the target for 
synthesizing SHN. In Fig. 14 we give evaporation residue excitation functions of the reactions 40 Ar, 
50 Ti, 54 Cr, 64 Ni+ 238 U in the 2n-5n channels. The results show that the 4n channel in the reaction 
40 Ar+ 238 U has the larger cross sections with 2.1 pb at an excitation energy 42 MeV. This reaction 
is being used to synthesize the superheavy nucleus Ds with HIRFL accelerator at Institute of 
Modern Physics in Lanzhou. The reactions 50 Ti, 54 Cr, 64 Ni+ 238 U lead to the cross section smaller 
than 0.1 pb. The isotopic trends based on the U isotopes are also investigated using the DNS model 
as shown in Fig. 15. Calculations show that the isotopes 235 U and 238 U are favorable in producing 
SHN. The cross sections are reduced with increasing the mass numbers of the projectiles. Other 
reaction mechanisms to synthesize SHN have to be investigated with theoretical models, such as 
the massive transfer reactions, and the complete fusion reactions induced by weakly bound nuclei. 
Work in these directions is in progress within the framework of the DNS model. 

4 Conclusions 

Using the DNS model, we systematically investigated the production of superheavy residues in 
fusion-evaporation reactions, in which the nucleon transfer leading to the formation of the su- 
perheavy compound nucleus is described with a set of microscopically derived master equations 
that are solved numerically and include the quasi-fission of the DNS and the fission of the heavy 
fragments. The fusion dynamics and the evaporation residue excitation functions in the 48 Ca 
fusion reactions are systematically investigated. The calculated results are in good agreement 
with the available experimental data within the error bars. Isotopic trends in the production of 
superheavy elements are analyzed. It is shown that the isotopes 235 ' 238 U, 237 Np, 241 ' 244 Pu, 242 Am 
and 2 45,247,250 Cm5 248,249 Bk and 25i,252 Cf in the 3n cha nnels, and 230 Th, 244 Pu, 24 8> 2 50 Cm in the 4n 

channels are favorable for producing the superheavy elements Z=110, 112 and 113-118, respec- 
tively. The evaporation residue excitation functions of the reactions 40 Ar, 50 Ti, 54 Cr, 64 Ni+ 238 U in 
the 2n-5n channels and the isotopic trends with 40 Ar, 48 Ca, 50 Ti, 54 Cr, 58 Fe and 64 Ni bombarding 
U isotopes are also studied. 
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Table 1: Parameters used in the calculation of the level density. 
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Figure 1: Calculated values of the level density parameters as a function of the atomic mass. 
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Figure 2: Calculated capture, fusion and survival probabilities as functions of the relative angular 
momenta in the reaction 48 Ca+ 208 Pb at excitation energies of the compound nucleus of 20 MeV 
and 40 MeV, respectively. 
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Figure 3: Comparison of the calculated fusion-fission excitation functions and the available exper- 
imental data for the reactions 16 O+ 208 Pb, 16 0+ 238 U, 36 Ar+ 148 Sm and 26 Mg+ 238 U. 
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Figure 4: Calculated quasi-fission mass yields for the reactions 48 Ca+ 244 Pu and 48 Ca+ 248 Cm at 
excitation energies of the compound nuclei 42 MeV and 33 MeV, respectively, and compared them 
with the available experimental data [39J. 
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Figure 5: The calculated evaporation residue excitation functions with 232 Th, 238 U, 237 Np and 
243 Am targets in 48 Ca induced reactions, and compared with the available experimental data 
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Figure 6: The same as in Fig. 5, but for the targets 242 > 244 Pu and ' 248 Cm to produce superheavy 
elements Z=114 and 116. 
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Figure 7: The same as in Fig.5, but for the targets 247 Bk, 249 Cf, 254 Es and 257 Fm to synthesize 
superheavy elements Z=117-120. 
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Figure 8: Maximal production cross sections of superheavy elements Z=102-120 in cold fusion 
reactions based on 208 Pb and 209 Bi targets with projectile nuclei 48 Ca, 50 Ti, 54 Cr, 58 Fe, 64 Ni, 
70 Zn, 76 Ge, 82 Se, 86 Kr and 88 Sr, in 48 Ca induced reactions with actinide targets by evaporating 3 
neutrons, in comparison with available experimental data [H [2j HI Wl\ . 




Figure 9: Isotopic dependence of the calculated maximal production cross sections in the 3n 
evaporation channel and the corresponding excitation energies in the synthesis of superheavy 
elements Z=110 and 112 for the reactions 48 Ca+^Th and 48 Ca+' A U, and compared with the 
experimental data [3D], W2\ . 
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Figure 10: The same as in Fig.9, but for isotopic targets Np and Pu to produce superheavy 
elements Z=113 and 114. 
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Figure 11: The same as in Fig.9, but for isotopic targets Am and Cm to synthesize superheavy 
elements Z=115 and 116 in 3n and 4n channels. 




Figure 12: The same as in Fig.9, but for isotopes Bk and Cf in Ca induced reactions. 
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Figure 13: (a) the inner fusion barrier, (b) the fission barrier of the compound nucleus and (c) the 
neutron separation energy as a function of the mass numbers of the isotopic targets Cm in the 
reactions 48 Ca+ A Cm. 
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Figure 14: The evaporation residue excitation functions in the reactions 40 Ar, 50 Ti, 54 Cr, 

64 Ni+ 238TJ. 
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Figure 15: The production cross sections in the 3n channels as a function of the mass number of 
the isotopic targets U with projectiles 40 Ar, 48 Ca, 50 Ti, 54 Cr, 58 Fe and 64 Ni. 
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